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ABSTRACT 



U 

<^ . We show that the process of photoionizing a gas of atomic hydrogen 

and helium by line radiation whose energy is slightly above the helium 
single-ionization threshold is unstable if the helium fraction by number is less 
than approximately one half. However, in the two scenarios we consider here, 
based on the Decaying Dark Matter (DDM) model of cosmological reionization, 



■ there is no significant growth. In the first scenario we consider ionization and 

recombination to be approximately in equilibrium. This is relevant to high 
photon flux rates and early reionization, but in that case the heating is balanced 
by Compton cooling, which is very stabilizing. In the second scenario we ignore 
recombination. This is relevant to low photon flux rates or to the last stage of 



the reionization. In that case there is too little growth on a cosmological time 
scale to be significant. 



Subject headings: cosmology: dark matter — cosmology: early universe — 
cosmology: large-scale structure of universe — hydrodynamics — instabilities 



1. Introduction 



The hydrogen in the intergalactic medium (IGM) is highly ionized and has been at 
least since z = 4.3 flSteidel fc Sargent 1987| ; piallongo et al. 1994| ). Recent observations 



indicate that the helium in the IGM is mostly at least singly ionized, at least at large 
redshifts ( [Reimers fc Vogel 1993| ; |Jakobsen et al. 1994] ; |Miralda-Kscude 1993| ). Though the 
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IGM might be ionized mainly by quasars ( [Mciksin fc Madan 1993 ; IVladau 1992| ) or an early 



generation of stars ( [Fukugita fc Kawasaki 1994j ), additional ionization by photons emitted 
by decaying dark matter (DDM) has many appealing features ( jl'lannery fc Press 1 9 79| ; 
[Sciama 1990 ; [Sciama 1993| [Sciama, Kees, fc Scott 19UT| ). 



In particular, if a generation of neutrinos has a mass sufficient to close the universe, 
and if such a particle decays into a much lighter particle and a photon, the photon's energy 
will be near the ionization energies of hydrogen and helium. Specifically, 

£ 7 « -m v = ^91.5/i 2 eV (1) 

where h is the Hubble constant in units of 100 km/s/Mpc. Thus, for h > .55, the decay 
photons would have sufficient energy to ionize hydrogen and for h > .74, they could (singly) 
ionize helium. 

The decay photons would also provide a source of out-of-equilibrium energy which 
might lead to small-scale structure formation. [Hogan (1992|) proposed a mechanism of 



small-scale structure formation through the ionization of a pure-hydrogen plasma, but 
Bradford fc Hogan (1994|) showed that this mechanism does not work. We have continued 



to seek possible instabilities of the system, since a neutrino-dominated universe requires a 
non-gravitational mechanism for spawning structure on galactic scale and below and there 
is no fundamental reason why this system should be stable. The free energy per atom is 
sufficient for many e- foldings of instability growth (E 1 ^> kT), and the scale of instability 
can be comparable to proto-galactic structure. A hydrodynamical instability would mimic 
many of the desirable features found in the gas component of a Cold Dark Matter (CDM) 
model, as described for example in Miralda-Escudc fc Rccs (1993|). 



1.1. Helium-Ionization Instability Mechanism 

Meiksin (199^) suggested that there is an instability if DDM photons have sufficient 
energy to ionize helium. The instability is thermal, i.e. it is based on low-density regions 
being preferentially heated. Low-density regions receive more heat because they have a 
larger ratio of neutral hydrogen to neutral helium and a hydrogen ionization deposits 
more energy in the plasma than a helium ionization (because there is a larger difference 
between the photon energy and the ionization threshold). The ratio of neutral hydrogen to 
neutral helium is larger in low density regions because they are more ionized and the higher 
ionization cross-section of helium (for photons just above its ionization threshold) ensures 
that it will absorb photons out of proportion to its abundance. 
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The reason that low-density regions are more ionized is different in the two scenarios 
considered below. When recombination equilibrium holds, it is due to the fact that 
recombination happens in proportion to the square of the density (because a positive ion 
and an electron need to find each other) while the ionization rate is proportional to the 
density (the rate of production of photons is assumed to be unrelated to the density). When 
recombination equilibrium doesn't hold, the low density regions are more ionized because 
there are more photons per atom. 

Though any source of uniformly distributed photons will do for this instability, we 
concentrate on decaying neutrinos and parameterize the photon production rate by the 
neutrino lifetime in units of 10 24 seconds, r 24 . 



2. Perturbation Expansion and Matter Evolution Equations 

We expand to first order in small perturbations about background values, and expand 
the perturbations in a Fourier series. We normalize some perturbations by the unperturbed 
value, e.g. for the number density, n, 



n = n 



(l + 51nne <& *) (2) 



The zero superscript denotes the unperturbed variable; it will be dropped in subsequent 
formulas for simplicity. 

The evolution equation for the internal energy density, u, is 

Ft u = lu Ft lnn + v (3) 

where n is the number density, 7 is the adiabatic expansion coefficient (5/3 for a monatomic 
gas), and V is the heating and cooling term.f] This equation linearizes to 

d d 

VS In u + u—5 In u = •yu—6 In n + 5V + V u 6 In u (4) 
dt at 

The relationship between density and pressure (or internal energy) perturbations is 
(see ( Peebles 19801) ) 

——5\nn + 2H—5\nn = — -— 5\nu (5) 

dt dt dt a 2 7 w 



3 Table 1 lists the symbols used in this article. 
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where k is the comoving wave number, H is the Hubble constant, a is the cosmic scale 
factor and is equal to 1/(1 + z), and v s is the speed of sound (remember that the pressure 
is just 7 — 1 times the internal energy density). 

Here 

vj_ = po = (l + x)T = (7 - 1> (6) 
7 Po rn proton (1 + 3 a) nm proton (l + 3//) 
where /x is the fraction by number of helium nuclei and \ ls the ratio of free electrons to 



total hydrogen and helium nuclei. Note that m p f ° t °^^ — tl j s the average mass of a plasma 



™.proton(l + 3/j) 

particle (since m e i ect ron < m pro tcm)- 
Combining eqs. 01 and R gives 



d d d ^ 9 9 n 9 t> 2 fc 2 r ^ 

SSS sln " + a2 SS flE " + a 'S iln " + ^ p = (7) 
o, = 2h(-- p u ) + h 2 + v 2 X 



S o V / 

v u J cr 

a 2 =--V u + 4H (9) 

Since a x and a 2 are positive (note that P u is zero for ionization heating and is negative 
for Compton cooling), we will have an instability if j^— < 0, which we will show to hold 
for the ionization instability considered alone, though it doesn't hold when the ionization 
heating is balanced by Compton cooling. 

Since a\ and a 2 are positive, they lead only to damping of growth. We ignore them 
from now on, which can only be generous to the instability. 

The analysis to this point has been general to any sort of heating. The form of 5 In V 
will be different for the two instability scenarios here, though in both V depends directly 
only on \ an d the relationship between n and x is used to get the variation of V with n. 
For the "recombination equilibrium" scenario, we will find SV in terms of 8 Inn and the 
equation will be third order; for the "no recombination" scenario we will have SV in terms 
of / S In ndt and the equation will be fourth order. 



3. Scenario in which Recombination and Cooling are Important 

3.1. Photon Energy Levels 



At the heart of the ionization instability mechanism is the fact that a helium ionization 
deposits more energy than a hydrogen ionization. However, a significant part of that 
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difference is lost when the helium atom recombines, emitting photons which might then 
ionize a hydrogen atom. 

A direct recombination to ground will yield a photon at the ionization energy of He, 
which we label E\\ if this photon ionizes a hydrogen atom, no energy will have been lost in 
the helium ionization. A recombination to an excited state of helium will emit a low-energy 
photon which can't ionize hydrogen, so its energy is lost. All such recombinations cascade 
(emitting other low-energy photons) to one of three states, which then decays directly to 
ground, usually emitting a photon which can ionize hydrogen (so we need to keep track of 
it). The 2 3 P triplet state and the 2 l P singlet decay by emitting one photon at energy Ei or 
E 3 respectively. Helium atoms which end up in the He 2 X S state can only decay to ground 
by the emission of two photons. We label the 2 X S state energy E4 and the average energy 
of the two-photon-decay photons which can ionize hydrogen E$ (see eq. ( |A2| ) ) . Finally, we 
label the the primary photon energy (the energy of the decay photons) Eq. 

On average p\ = 0.554 hydrogen ionizations result from each two-photon decay of the 
2 1 jS i state (eq. ( |A1Q ); all other helium recombinations lead (in the optically thick limit) to a 
hydrogen ionization. 



See Osterbrock (198S) for more background on these radiative processes. 



3.2. Ionization and Recombination Equations 

The heat energy absorbed will be proportional to the numbers of photons in each level 
times the absorption cross section for those photons: 

V = n 7 o[(£'o - Ry)cr H o"-Hi + (E - -Ei)a He o^Hei] + n 7 i(#i - Ry)er H in H i + 

n j2 (E 2 - Ry)a m n m + ^73(^3 - %Vh3™hi + n j5 (E 5 - Ry)a m n m . (10) 

^hi, tt-hiI) ^Hei; n Heii are the number densities of neutral and singly ionized hydrogen and 
helium atoms. n 70 , n 7 i, etc., are the number densities of photons with energy E , Ei, etc. 
Ry is one Rydberg, the hydrogen ionization energy. The sigmas are the absorption cross 
sections for the given substances with photons of the indicated energy, e.g. o"Hei is the cross 
section for the ionization of helium by photons with energy E\. 

The following equations express the assumption that the creation and destruction of 
photons at each level are in equilibrium: 

= (o"HeO^HeI + O"H0^Hl)^ 7 
ttHel^Hell^e = (^Hel^He I + 0"H1^H l)« 7 l 
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ttHe2 3pn H eII^e = 0"H2^HI^ 7 2 
O^^HeH^e = 0"H3^HI^ 7 3 
aHe2i5^HeII^ePl = 0"H5™HI™ 7 5 (H) 

9ft is the rate of production of primary photons per unit volume (in the decaying neutrino 
model it is n v /r v ). The a's are the net recombination coefficients into various states, as 
defined in Table 1. Direct recombination of hydrogen to ground is ignored since (in the 
optically thick limit) the photon produced immediately ionizes another neutral hydrogen 
atom. 

The following equations say that the creation and destruction of each species of ion are 
in equilibrium: 

amnnune = (cr H o™ 7 o + 0"Hin 7 i + cr m n l2 + <?"H3™ 7 3 + o- m n^)n m 

"He^Hel^e = (o"HeO™ 7 + 0"Hcl™ 7 l)™HeI (12) 



We also use the following relationships (from Ostcrbrock (1989)) ) between recombination 
coefficients: 



«He = 


«Hel + «He2 




3 


«He2 3 P = 


^«He2 






C*He2 IP = 








«Hc2 IS = 


l (r He2 ) 



(13) 



3.3. Analysis 

Combining eqs. (|ll|), fll2"D and (|H|) gives the density and hydrogen ionization coefficient 
in terms of the helium ionization coefficient: 

(1 - /i)cr H lQ!He + hiPl + H)0- He lA*(l - XHe)«He2 ,x 

XH = J-, w J-, v ; v XHe (,14) 

(1 - /i)(CT H cl(l - XHe)aH2 + CTHlXHettHe) 

n 2 = ggelg(l - XHc) , 15 x 

((1 - fl)<J m (l - X H )«He + 0"HelA*(l ~ XHe)aHe2) ((1 ~ A*)XH + HXKe)XRc 

Applying those relations to the formula for V (eq. (|T(J|)) gives 

v = k(ae a «f a ™ . eJ] \^4r^ ( 16 ) 

V CTHel/^aHe2 + 0"Hl(l ~ fl)a He 7 1- D^He 
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^ (AJ5 - -E L oss)0-HeiaH2 + (-^Loss ~ ^(1 ~ Pl)ASj CT H1 a H e « H e2 ^ 

(AE - £ , Loss )cr H ei^aHe2 + AEcr H i(l - /i)a H e "H2 

CT H eiaH2 - 4(1 -Pl)0"HlttHe « H e2 /.. D x 

^ = n ^ 1 ^ 

0"Hl(l - /iJttHe + C r HelAtaHe2 «H2 

(19) 

A£" is the difference between Eq and a Rydberg and -Eloss is the difference in the energy 
deposited when a primary photon ionizes a hydrogen atom instead of a helium atom: 

£loss = E x - ^(E 5Pl + 2E 3 + 9E 2 + Ry(l - Pl )) = 4.97 eV. (20) 

To get positive feedback in the heating, we need the heating to go down with increasing 
density and hence to go up with increasing ionization. For that we need A < B or 

CT He l(l - /i)«H2 > 0"Hl(l - |")aHe + T^(pi + ll)<7 H eljUO!He2 (21) 

which holds since /i ^ ~ (i.e. there is more hydrogen than helium), am <C CHei (i-e. the 
absorption cross-section for helium is larger at this energy), and |aH e ~ «He2 ~ «H2 (i-e. all 
the recombination coefficients are comparable). 

Filling in the numeric values given in Table 1 and taking E = E\ (that is, ionization 
photons right at the helium ionization threshold, which is most generous to the instability), 
we get V and n 2 in terms of the variable XHe and the parameter T24: 

eVz' 3 1 - 0.180xHe 
cm 3 s 1 - 0.246xhc 

ot) x ur*— (1 - 0.736 X He) 2 , . 

r 24 cm 6 (1 - 0.243 X Hc) (1 - 0.326 X Hc) XHe 2 



V = 1.12 XlO- 21 ^ ^- U -^He (22) 



where 2;' = 1 + z. When varied these give 

5V = 0.0662 V- r; -5v He (24) 

(1-0.180xhc)(1 - 0.246xhc) A 1 J 

(note the small coefficient in this equation, indicating the weakness of this instability) and 

5 = Sn (1 - 0.246x He ) 2 (l ~ 0.351 X He) 2 (l - 0.736xhc) 2 Xhc 

XHc n (1 - 0.600xHe + 0.266xhc 2 ) (1 - 0.243 X Hc) (1 - 0.295 X h c ) (1 - 0.326xh c ) (1 - 0.771 X He) 

(25) 

Defining the coefficient multiplying 5 In n in eq. [7] as 

v 2 k 2 

a = ^r~V n (26) 

GT7W 
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and using a 2 = (H/H ) 3 we get 

2.56(1 - 0.35xhc) 2 (1 - 0.74xHe) 2 XHe(Mpcfc) 2 # 3 



a 



1 - 0.60XHe + 0.23XHc 2 )(l - 0.24 X He)(l - 0.29 X He)(l - 0.33 X He)(l - 0.77 X Re)z' 5/2 T 24 

(27) 



We can use eq. ^ and the known background value for the density, n = 1.04 x 10 7 
(this assumes r\ = 3 x 10 10 ) to get the background ionization for a given redshift: 

r MZ ' 3 = 4.06x10' £ ~ °y ^ (28) 

XHe 2 (l - 0.27XHe)(l ~ 0.37x H e) 

Using this equation to eliminate the redshift gives 

9.30 x 10~ 5 (1 - 0.35 X Hc) 2 (l - 0.74xHc) 1/3 XH / e 3 (MpcA;) 2 i7 3 



a 



(1 - 0.60 X He + 0.23 X He 2 )(l - 0.24 X Hc) 1/6 (1 - 0.29 X Hc)(l - 0.33xhc) 1/6 (1 - 0.77 X He)r 2 1 4 /6 

(29) 

This is largest (in magnitude) at xh = 1, where it is 



a = -2.79 x 10~ 4 V P J 6 (30) 



21 



For a sufficiently small time, a will be effectively constant and the growth rate will be 

Ug = (31) 

We will show that there is no growth for any such small time interval, so there will be no 
growth for any larger time interval. 



3.4. Range of the Instability 

The minimum scale for this instability is the point where the assumption that the 
medium is optically thick breaks down; it is the optical path length (which should be the 
minimum of the lengths for a helium or a hydrogen ionization, but for convenience we take 
it to be the length for a helium ionization): 

^min 7Z \ (32) 

CT H e(l - XHeJ^He 



or (still with E = E% 



l min z' = 5.46Mpc- —2 (33) 

I 1 — Xn e )z 
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or in terms of wave number, 

U = t^7 = l.lSMpc-^l - xu c )z' 2 . (34) 

The maximum scale comes from requiring that the growth rate be larger than the 



expansion rate: uo„ > H. Taking ui g from eq. 31 gives 



k min = 104. Mpc -1 -^ 1 / 12 x 



(1 - 0.60XHC + 0.27 X He 2 )(l ~ 0.24 X He) 1/6 ( 1 ~ 0.30XHe)(l ~ 0.33x H e) 1/6 (1 ~ ^^X^ 



\ (1 - 0.35 X He) 2 (l - 0.74 XH e) 1/3 XH / e 3 



where the ionization is related to the redshift through eq. . 

Requiring that the instability have a non-zero range of scales, or 

^min *^ k mSuX (36) 

gives a limit on xne or equivalently z' for a particular r 2 4. Figure |I| shows the lower bound 
on z' as a function of T24. The lowest z' at which there is a finite range for the instability is 
z' = 21.0 (for r 24 = 20.7 and XHe = 0.665). 

An additional constraint is that the assumption of ionization equilibrium be valid. 
Define the ionization rate per helium ion as u;R ec = an e n e = a^eXP- F° r these parameters 
a>Rec = 1.01.ff so equilibrium still (marginally) holds. 



3.5. Compton Cooling 

In the post-recombination, pre-galactic universe of our model, the dominant cooling 
mechanism is Compton scattering off the background radiation (since in the DDM model 
there will always be a significant ionized fraction). 

The power lost to Compton cooling is 

£ = AaT f Th ° mPSOn n e (T - T 7 ) = 4flT 4 Thompson _ ^ } ^ 

which consists of a cooling term proportional to «^ and a heating term proportional to n%. 
The cooling term increases with ionization and thus opposes any instability. The heating 
term is approximately constant with density or ionization (since x varies approximately as 
1/n) but is also (weakly) stabilizing. 
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Since Compton cooling depends more directly on the ionization (C x ~ C), its stabilizing 
effect dominates if the cooling is comparable in magnitude to the photo-ionization heating. 
With a simple model for the temperature^], we calculate that at the lowest redshift possible, 
z' = 21.0, the ratio of the cooling density dependence to the heating density dependence is 

r 

< -2.68. (38) 

' n 

Figure shows ^ as a function of T24 at the minimum z' using the same thermal model. 

1 — 1 ' n 

It's smallest magnitude is 

r 

< -2.00 (39) 

at = 24.3, r 24 = 55.7, and xtte = 0.409.Q 

Thus, the stabilization of Compton cooling will always dominate the destabilizing 
effects of the heating. 



4. Scenario in which Recombination is Ignored 

At redshifts or ionization rates which are low enough that Compton cooling is not 
significant, recombination is also less important. We model this regime as having no 
recombination, which strengthens the instability, but reduces the time over which it can 
operate. 



4.1. Ionization Equations 

We solve for the ionization by still assuming photon production to be in equilibrium 
with ionization: 

3? = (tfHeO^Hel + CT m n m )n j0 (40) 

but ionization is not in equilibrium with recombination: 

dXUe 



dt 



n l0 a Hc0 (l - XHe) (41) 



4 To guarantee that the model temperature was a lower bound, we took the ionization in C to be constant 
at 1 and the ionization in V to be constant at and we approximated the resulting integral for u in a way that 
gives about half the actual value for the difference in the matter temperature and the photon temperature. 

For these parameters, wr cc = Q.Q12H , so the assumption of ionization equilibrium is invalid. Enforcing 
that assumption would require an even less favorable ratio. 
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cit 

Taking the ratio of eqs. (O) and ( |4"2"D we get 

dxn cho 1 - Xh 



'' v " n 70 (THo(l - Xh)- (42) 



(43) 



<^XHe CTHeO 1 — XHc 

or 

1 - Xh = ci(l - Xhc) CTHc0 (44) 



where ci is determined by the initial conditions. This is an exact equation which holds 
for the background ionization and which can be varied to yield a relationship between the 
hydrogen and helium perturbations: 



<5xhc = ci— ^(1 - Xhc) CTHc0 ^Xh (45) 

°"HeO 



From eqs. (gTJ), ([PJ), and (0) we get 



3? 

;i-/i)(i-x H ) + Mi-XHe) = - / (46) 



n 



Note that - is a constant in the unperturbed background, so for background quantities: 

(l-/i)(l- XH )+Ml-XHc) = 1 (47) 

n 

which gives the ionization as an implicit function of time or (as we shall use it) gives the 
time as a function of the background ionization. 

Using the following formula (valid for a flat, matter-dominated universe) we can express 
the redshift in terms of time and thus ionization. 

t = t^r 1 (48) 

6 tf 



Varying eq. EBI gives 

3? 



(1 - fi)6xn + fiSxuo = — J Slnndt (49) 
which relates the variation of the density to the variation of the ionization. 
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4.2. The Variation of Heating with Density 



The heating function, V, is simpler than in the first scenario because there are no 
photons at reprocessed energy levels: 



V = n j0 [(E - Ry)cr m n lu + (E - E^ciHeo^Hei]- 



(50) 



Using eq. [40] to eliminate n 70 gives 



v = n 



(E - Ry) - (E 1 - Ry) 



I I l-M °H0 1~XH 
A» °"Ho() 1— XHe 



(51) 



As in the equilibrium case, the heating function depends on the density only through 
the ionization. Varying V with respect to the ionization and using the relationship between 
ionization and density variations (eq. £49]) gives 5V n which substituted into eq. [7] gives 



Odd a a a f 

T^^T^$\nn + a 2 T^T^51nn + a 1 —5lnn-a / Shindt 
at at at at at at J 



d d 



d 



(52) 



where 



(7- cr H0 / a m \E 1 -Ry 1 - Xn 

0-0 — — ~ I 1 



1 



ft' 



(3/i + l)yU 2 a He V~ (ThcO/ "^proton (1 ~ XHe) 2 ( \ + l^M^HO. 1~Xh ) 3 U 

and a\ and a<i are given by eqs. § and P. 



2^/2 



k z' 



(53) 



4.3. Bounds to the Growth 

In this scenario an instability can grow, though the total growth is less than one 
e-folding so it is insignificant. 

The time available for growth is bounded by the neutral atoms being used up, since 
there is no recombination. The starting ionization will be greater than zero since some 
ionization will occur while Compton cooling damps any growth. We take Xorig to be the 
recombination-equilibrium helium ionization (eqs. [15] and [14]) at the red-shift where the 
damping from the Compton cooling just becomes weaker than the instability growth. This 
is generous to the model in three ways: 1) assuming recombination equilibrium always over 
estimates the remaining neutral fraction; 2) the Compton cooling is calculated for matter 
at the photon background temperature, where it is lowest (the net cooling will be zero at 
this temperature, but there is a non-zero derivative of the cooling with respect to helium 
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ionization); and 3) the growth will actually be zero at that point, not its full value without 
Compton cooling. Figure |3] shows Xong as a function of r 2 4. z' Qrig and the original xh can 
be obtained from Xorig and t 2 4 using eqs. 128] and O. Those parameters then determine C\ 



from eq. 03 and the constant of integration in eq. 



The largest ionization, XFinai is the value of the helium ionization at which the 
mean-free-path exceeds the scale under consideration (i.e. it depends on k). It is determined 
for a given k from eq. [34] (given the redshift /ionization relationship of eq. 09]). 



To determine the total growth, we broke the evolution down in steps with ao taken to 
be constant at its maximum value in each step (the maximum is achieved for the smaller 
ionization at the start of the step). For constant ao (and ignoring a 2 and 03, which only 
damp any growth) the growing solution is <51nn cx exp( 1 */aot). We took steps in xne instead 
of time since we have time and redshift as a function of XHe and not vice-versa. Thus the 



growth in 5 Inn is bounded by exp(J2iL 1 (t(xi+i) — t(Xi)) y a o(Xi) f° r any N (with a tighter 
bound for a larger N) with Xi = Xorig + ^(XFinai - Xorig)- 

Figure [| shows the log of the total growth (for N = 40) as a function of Xorig and 
XFinai- It's largest value is exp(.83) at Xorig = -44 and XFinai = -79 (for these parameters the 
model is actually invalid because recombination is large, but it shows that within the space 
of validity the growth will be even less). 



5. Conclusions 

Though there is an interesting instability involving helium ionization by line radiation, 
it doesn't seem to have a cosmological application in the DDM scenario. As other candidate 
radiation sources presume the existence of small-scale structure, we conclude that this class 
of instability is unlikely to play a role in the initial formation of structure from smooth 
cosmic gas. 
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A. Calculation of Two-Photon Decay Parameters 

The He 2 X S state decays by the emission of two photons whose energy sums to E±. 
The strength of the transition which emits one of the photons in the range v to v + dv is 
A{y)dv. We did an empirical fit to the average of the two tables of values for A given in 
IJacobs (1971|) . Figure [5] shows the fit and compares it to Jacobs's values. 



The total strength for emission of two photons is A^ ot = J^ 4, A{y)dv. The total 
strength for the events in which a particular one of the photons is energetic enough to 
ionize hydrogen is A^y = A{y)dv. Since we have a hydrogen ionization if either photon 
is energetic enough, the proportion of decays which ionize hydrogen isQ 

Pl = 2—^ = 0.554 (Al) 

^Tot 

Weighting the strength of decay by the energy of the decay gives the average energy of 
emission, here taken over only those photons which can ionize hydrogen: 

Af v 4 vA(v)dv 

E 5 = JRy V ' = 16.1eV (A2) 
^Ry 

Since the absorption cross section depends on frequency, it should also be calculated as 
a weighted average: 

^hs = — - A (A3) 

A Ry 

= im0a H (E 5 ) (A4) 

so the weighting gives a 5% correction compared to using the absorption cross section for 
line radiation at energy E 5 . 



B. Ignored Effects 

Recombinational cooling, or the thermal energy lost when two particles which 
both have average kinetic energy T combine to one particle with average energy T, is 
insignificant. Since each ionization results in a recombination in the first scenario, the ratio 



6 This differs siightly from .56 given in Ostcrbrock (1989| ), perhaps because his value was not updated 
when he updated £4 from 20.7 to 20.6. 
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of recombination cooling to ionization heating is roughly the ratio of the average energy 
deposited by an ionization to the temperature, or about T/6.5eV. Compton cooling keeps 
the matter temperature comparable to the radiation temperature, which is much less than 1 
eV. For the second scenario there is no recombination, and thus no recombinational cooling. 
Note that no collisional ionizations happen at the low temperatures considered here. 

Collisional Transitions between helium levels are insignificant at the cosmological 
densities which are of primary concern in this analysis. The critical density given in 
Osterbrock ( 1989]) above which one must consider collisional transitions is ~ 5000cm -3 , 
much larger than cosmological densities of ~ 10~ 7 z' 3 cm~ 3 < 100cm -3 for z' < 1000. 

The magnitude of the photon redshift over one mean- free-path (for, e.g., an Ei 
photon to ionize helium) is 

Az = 1.8 x 10 -3 /i — a — (Bl) 

Z' - XHe) 

An Ei photon will be redshifted by more than one helium thermal Doppler over one 
mean-free-path for z ^ 100; this effectively reduces the absorption cross section of helium, 
which weakens the instability. The redshift over the mean-free-path of a one Rydberg 
photon is greater than a hydrogen thermal Doppler width for all z <; 50; for smaller redshifts 
than that, some hydrogen recombinations direct to ground will be allowed, which slightly 
strengthens the instability, as can be seen by considering eq. (pT[) with increased a^2- 



A distribution of initial photon energies will result in reduced growth over a line 
spectrum with an energy just above the Helium ionization threshold simply because any 
photons much above the threshold are less effective and any photons below the threshold 
actually oppose the instability. 

Lyman-alpha trapping will occur, since for z ^ 40, the redshifting in one mean-free- 
path is insufficient to bring a Lyman-a photon more than one thermal Doppler width away 
from the line center, but because the time photons spend traveling between encounters is 
much larger than the time a hydrogen atom spends in the 2 P state, only approximately 
10~ 7 of the neutral hydrogen atoms are excited, and this does not affect the instability. 

Photon diffusion (relaxing the assumption that the medium is optically thick) 
opposes an instability, since the more-ionized/less-dense regions are more optically thin and 
would thus lose photons (and therefor energy) to the less- ionized/more-dense regions. 

Thermal conduction opposes any instability and its effects are insignificant on the 
large scales considered here. 

Collisional ionization goes up with temperature, but it goes up with density squared 
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so it would also act against an instability. Its rate is insignificant at the temperatures and 
densities of the post-recombination universe. 

The temperature dependence of the recombination coefficients does not have a 
large effect because only the ratios of recombination coefficients enter the formulas, and the 
recombination coefficients for helium and hydrogen both depend on temperature to roughly 
the same power (3/4 for T « lOOOOfT). 

Any delay in thermalization has no effect, since Compton cooling and pressure are 
both linear in the sum of the energy of the particles, so the distribution of the energy does 
not matter. 
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Fig. 1. — Minimum z' = z + 1 vs T24 = t u /10 24 . For small T24 (large photon production 
rate), the minimum redshift is large because the large equilibrium ionization leads to a large 
mean free path. For large T24 the minimum redshift is large because the low ionization leads 
to low instability growth. 




20 40 60 80 



Fig. 2. — ^ vs T24. The curve follows the inverse of the minimum z' curve (since Compton 
cooling is much more effective at large z), with a shift towards larger r 2 4 (where the ionization 
and thus the cooling is less). 
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Fig. 3. — Starting xtie f° r no- recombination conditions vs r 24 . This is the recombination- 
equilibrium ionization where the Compton cooling ceases to dominate the heating. The steep 
slope near r 2 4 = reflects the fact that the heating approaches zero as xne goes to one, so 
there needs to be a high photon production rate to compensate. 
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Fig. 4. — Logarithm of the total growth factor vs the initial and final ionization for no- 
recombination growth. Since ionization increases with time, we clearly need the final 
ionization to be larger than the initial ionization, with larger values of their difference giving 
more time for growth. However, very small values for the original ionization correspond to 
low photon production rates, which give small overall growth. Large values for the final 
ionization correspond to large spatial scales, where the growth is also reduced. 
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Fig. 5. — He 2 1 S' two-photon decay rate versus energy (x = v/E±). The fitting formula is 
A(x) = -0.98p{2,x) + 2.83p{3,x)-A.26p{4,x) + 2.15p{5,x), where p(j,x) = 300((2a; - l) 2 ^' - 
1 — j((2x — l) 2 — 1)) (this form is chosen so that each term is symmetric around x=0.5 and 
is zero and has zero slope at x = and x — 1). The RMS error of the fit is 0.34. 
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Table 1. Symbols 



Symbol/ Value Explanation 

h The Hubble Parameter in units of 100 km/s/Mpc. 

Ug = —iu Instability growth rate. 

k Wave vector in Fourier decomposition. 

k = \k\ Wave number (inverse of size scale). 

n Total spatial density of hydrogen plus helium nuclei 

(neutral and ionized). 

u = ^rr(l + X)nT m Internal thermal energy of the gas. 

7 = 5/3 Adiabatic expansion coefficient. 

V Heat input per unit volume. 

v s The unperturbed speed of sound. 

\i = 4 _£ Y pa 0.073 Fraction by number of helium atoms (Y p m .24 is 

used here). 

w u = V u Miscelaeous rates. 
u k = v s k 
cj n = 

Pi = 0.554 Number of photons produced in a two-photon decay 

of the He 2 X S state which can ionize H. 

n 7 o,n 7l , etc. Spatial density of photons with energy E ,Ei, etc. 

^Heii — nnxue Number density of ionized helium. 

^Hei — nfi(l — XHe) Number density of neutral helium. 

nun = n(l — fi)xn Number density of ionized hydrogen. 

nm = n(l — /i)(l — xh) Number density of neutral hydrogen. 

n e = ^Heii + ^Hii = nx Number density of electrons. 

XHe Fractional ionization of helium. 

Xh Fractional ionization of hydrogen. 

X — f^Xiie + (1 — aOxh Fractional ionization of all species combined. 

<Th(E) = 6.30xl0~ 18 x Absorption cross section of hydrogen for photons 

x i-e-^y — cm Wlth ener gy E 
X = Ry 

y = vhi 
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Table 1 — Continued 



Symbol/ Value Explanation 

°Hm = CH(-E'm) Absorption cross section of hydrogen for photons 

with energy E m , for some m. 
ffHem = &Hm x Absorption cross section of helium for photons with 

(6.53^ - 0.22) energy E m , for some m. c 

«H2 = 2.59 x lCT 13 ^- Recombination coeff. for H to all excited levels. a 
«He2 = 2.73 x 10~ 13 ^- Recombination coeff. for He to all excited levels. a 
QHc2ip = | Q«Hc2) Recombination coeff. for He to excited singlet levels 

which cascade to the 2 l P state. 
«He2 1 5 — I (|o!He2) Recombination coeff. for He to excited singlet levels 

which cascade to the 2 l S state. 
«He2 3 P = |ttHe2 Recombination coeff. for He to all excited triplet 

levels (all of which lead to the 2 3 P state). 

«Hei — 1-59 x 10 _13 ^2_ Recombination coeff. for He to ground. a 

3ft = n u /r u Rate of production of primary photons. 

5 = ^ — 1 Dimensionless excess energy of ionization photon. 

6 = ^ Ratio of the matter temperature to the radiation 

temperature. 

z' — 1 + z Inverse of the cosmic scale factor. 

T24 = T„/10 24 sec Neutrino decay rate in units of 10 24 seconds. 



a From Osterbrock 1989 for T = 
b From Spitzer 1978. 
c From Brown 1971. 



10000K. 
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Table 2. Energy Values 



Sym. 


Value 


Explanation 


Eq 




Primary photon energy 


E 1 


24.6 eV a 


Ionization energy of He 


E2 


19.8 eV a 


Energy of He 2 3 P 


E3 


21.2 eV a 


Energy of He 2 X P 


E 4 


20.6 eV a 


Energy of He 2 *S 


E 5 


16.1 eV b 


Avg. energy of 7 absorbed by H in 27 decay of He 2 VS 1 


Ry 


13.6 eV 


H ionization energy 


AE 


E — Ry 


Energy deposited in a direct H ionization 


AE Uc 


Eq — Ei 


Energy deposited in a He ionization 



a From Osterbrock 1989. 
b See Appendix C. 
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